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For topologically nontrivial and very narrow bands, Coulomb repulsion between electrons has been 
predicted to give rise to a spontaneous fractional quantum-Hall (FQH) state in absence of magnetic 
fields. Here we show that strongly correlated electrons in a t2 g -orbital system on a triangular 
lattice self-organize into a spin-chiral magnetic ordering pattern that induces precisely the required 
topologically nontrivial and flat bands. This behavior is very robust and does not rely on fine tuning. 
In order to go beyond mean field and to study the impact of longer-range interactions, we map the 
low-energy electronic states onto an effective one-band model. Exact diagonalization is then used 
to establish signatures of a spontaneous FQH state. 



The Integer Quantum Hall (IQH) effect pQ is a prime 
example of an electronic state that cannot be classified 
within the traditional framework of symmetry breaking, 
but is instead characterized by a topological invariant [2] . 
It is by now theoretically well established that an exter- 
nal magnetic field is in principle not needed and that 
states within the same topological class as IQH states 
can be realized in lattice models, if time-reversal sym- 
metry is broken by other mechanisms, e.g., by complex 
electron hoppings [S]- Related topologically nontrivial 
Quantum Spin-Hall (QSH) states even occur in systems 
where time-reversal symmetry is not broken at all [1HH]j 
see Refs. [9] [TUj for reviews. At present, many intrigu- 
ing features intrinsic to topologically non-trivial states 
have been observed in the absence of magnetic fields, 
such as the metallic Dirac cones at the surface of a topo- 
logical insulator [TTJ [12], or the QSH effect in quantum 
wells [H[T1]. 

Fractional Quantum Hall (FQH) states [15] are topo- 
logical states that can be seen as composed of quasi- 
particles carrying an exact fraction of the elementary 
electronic charge |16j . Apart from the fundamental in- 
terest in observing a quasi-particle that behaves in many 
ways like a fraction of an electron, some FQH states also 
have properties relevant to fault-tolerant quantum com- 
putation [T7]. Very recently [TM2"U] . it was suggested 
that lattice-FQH states may similarly arise without a 
magnetic field, in fractionally filled topologically nontriv- 
ial bands. 

In contrast to the IQH and QSH effects, which can be 
fully understood in terms of non-(or weakly-)interacting 
electrons, interactions are an essential requirement for 
FQH states, which places demanding restrictions on can- 
didate systems: One needs a topologically nontrivial 
band that must be nearly flat - similar to the highly 
degenerate Landau levels - so that the electron-electron 
interaction can at the same time be large compared to 
the band width and small compared to the gap sepa- 
rating it from other bands [T8Tf2"0"] . If the requirements 
can be fulfilled, however, the temperature scale of the 
FQH state is set by the energy scale of the interaction. 



This can allow temperatures considerably higher than 
the sub-Kelvin range of the conventional FQH effect, 
which would be extremely desirable in view of potential 
quantum-computing applications. Moreover, the lattice 
version of FQH states [21] may have unique and different 
properties. [22) . 

In most recently proposed model Hamiltonians [T51 - 
I2TJ1 [2SH2S] , the topological nature of the bands was intro- 
duced by hand and model parameters have to be care- 
fully tuned to obtain very flat bands. As potential real- 
izations, "purpose built" physical systems in oxide het- 
erostructures [25] or optical lattices [T!5] were suggested. 
On the other hand, topologically nontrivial bands can 
in principle emerge spontaneously in interacting electron 
systems [27] [28], e.g., for charge-ordered systems [29] [30] 
or for electrons coupling to spins in a non-coplanar mag- 
netic order [31] [32]. We demonstrate here that such a sce- 
nario indeed arises in a Hubbard model describing elec- 
trons with a t2 g orbital degree of freedom on a triangular 
lattice: a ground state with topologically nontrivial and 
nearly flat bands is stabilized by onsite Coulomb interac- 
tions. Upon doping the flat bands, longer-range Coulomb 
repulsion induces FQH states. 

t 2 g orbitals on the triangular lattice. — The build- 
ing blocks of our system are oxygen octahedra with a 
transition-metal (TM) ion in the center, the most com- 
mon building block in the large and versatile material 
class of TM oxides. Local cubic symmetry due to the 
oxygen ions splits the d-orbitals into t2 g and e g lev- 
els, and it has been shown that orbital degrees of free- 
dom of either kind can substantially reduce the width 
of topologically nontrivial bands 24J. Here, we concen- 
trate on the t2 g orbitals illustrated in Fig. [ija), which 
are further split by a crystal-field due to the overall 
lattice geometry. On a triangular lattice, we find one 
ai g and two e' g ± states, see Fig. (lib), with a splitting 
ifjx = A JT (n eg+ + n Eg _ - 2n aig j/3 depending on the 
Jahn- Teller effect and the lattice |33J. Electron hopping 
along nearest-neighbor (NN) bonds consists of terms t via 
ligand oxygens and tdd due to direct d-d overlap [331 EI] j 
hopping matrices are given in [35) . We set here n < 3 and 
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FIG. 1. (Color online) Triangular perovskite lattice and ti g 
orbitals. Oxygen octahedra are indicated by lines, with black 
lines illustrating the front facets. Thick dotted (dashed, solid) 
lines indicate nearest-neighbor bonds along lattice vector ai 
(a2, 83). (a) Shows two d xy orbitals (top) and one d xz and d yz 
orbital (bottom). In (b), the orbitals reflecting the three-fold 
lattice symmetry are shown: The two e' g orbitals (bottom), 
which differ by their complex phases, will turn out to be half 
filled, while the ai g orbital (pointing out of the plane, see top) 
forms nearly flat bands with non-trivial topological character 
that can support spontaneous FQH states. 
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FIG. 2. (Color online) Spin-chiral magnetic phase with topo- 
logically nontrivial bands stabilized by onsite Coulomb inter- 
actions in t zg electrons on a triangular lattice, (a) Chiral 
magnetic order, the sites of the unit cell are labeled by I to 4. 
(b) The spins on the four sites can be seen as pointing to the 
corners of a tetrahedron, i.e., the pattern is non-coplanar. (c) 
One-particle energies on a cylinder (periodic boundary con- 
ditions along x) in the mean-field [35] ground state of the 
t2 g multiorbital Hubbard model, which is given by the pat- 
tern shown in (a). States drawn in black (grey) have more 
(less) than 33% a lg character, dashed and dotted lines indi- 
cate edge states with more than 33% of their weight on the 
top (bottom) row of sites. The arrows f (\.) indicate states 
with electron spin mostly (anti-)parallel to the local quanti- 
zation axis, which can be seen as the lower (upper) Hubbard 
band. The filling is 2.5 electrons per site, slightly less than 
half filling. Parameters used were t — 1, tdd ~ 0, U/t = 12, 
J/t — 3, Ajt/£ = —6. The figure of merit M, which is given 
by the ratio of the gap separating the two a\ g subbands of 
the lower Hubbard band and the band width of the highest 
subband of the lower Hubbard band, is M ~ 14. 



choose t > [33] as unit of energy, but analogous results 
hold for n > 3, t < 0, and tdd —> —tdd, A JT — » — A JT due 
to particle-hole symmetry. 

In TM oxides, Coulomb interaction is substantial com- 
pared to the kinetic energy of t 2g orbitals and spin-orbital 
physics induced by correlations are known to be rich in 
tig systems on triangular lattices [3H IMj ■ We take into 



account the onsite interaction including Coulomb repul- 
sion U (intra-orbital) and U' (interorbital) as well as 
Hund's-rule coupling J. We employ a mean-field ap- 
proximation with a decoupling into expectation values 
of densities (ni !ajCr ) — (c\ aa .c iag .) for site i, orbital a, 
and spin a [37J 13E] • The spin is thus reduced to its z- 
componcnt m\_ a = (n.i,a,t — ni tOL ^)/2 and non-collincar 
magnetic patterns are treated by allowing for a site- 
dependent spin-quantization axis expressed by angles 0j 
and 4>i. The change in quantization axis from site to site 
manifests itself in a complex Berry phase for the hopping 
terms [39] . Numerical optimization is used to find the 
0i and 4>i giving the magnetic ground state, permitting 
arbitrary magnetic orderings with unit cells of up to four 
sites, including all phases considered in Ref. (30J. For 
simplicity, we present here results for J/U = 1/4 and the 
relation U' = U — 2 J between the Kanamori parameters 
was used, but we have verified that the results presented 
remain robust for other choices. For details see [35] , 

For wide parameter ranges (see below), the ground 
state is the non-coplanar spin-chiral phase illustrated 
in Fig. [2|a,b). As demonstrated in the context of the 
Kondo-lattice [38l [40] and the Hubbard [38l [41] mod- 
els, this magnetic order leads to topologically nontrivial 
bands, which can also be seen in the one-particle bands 
shown in Fig. [2jc) . The chemical potential lies within the 
a\ g states of the lower Hubbard band, where the electron 
spin is mostly parallel (labelled by f) to the direction 
defined by the spin-chiral pattern. Dashed and dotted 
lines decorate states living on the top and bottom edges 
of a cylinder, they cross the chiral gap exactly once as 
one left- and one right-moving edge mode, indicating the 
different Chern numbers associated with the two bands 
directly above and below the chemical potential. Such a 
spontaneous IQH state is already rather exotic and has 
recently been shown to support fractionalized excitations 
bound to vortices |4"2"] . 

Figure [2jc) also indicates that the upper chiral sub- 
band has a very small width, ~ 14 times smaller than 
the chiral gap. One can quantify the band flatness by 
a figure of merit M given by the ratio of the gap to the 
band width. Its dependence on various parameters of the 
Hamiltonian is shown in Fig. [3| It peaks at M > 40, but 
the more striking observation is that it is above 5 or even 
10 for wide ranges of U, Ajx and tdd, in contrast to many 
other proposals that require carefully fine-tuned parame- 
ters 12^25 1 115] . Nearly flat chiral bands are thus 
very robust in this system and both their topological 
character and their flat dispersion emerge spontaneously 
with purely onsite interaction and short-range hopping, 
without spin-orbit coupling or any explicit breaking of 
time-reversal symmetry. 

Mapping to an effective model. — For large onsite in- 
teractions and large crystal field splitting U, J, |A JT | >• 
t, tdd, the three-orbital model with fillings between 2 and 
3 electrons per site can be mapped onto the one-band 
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FIG. 3. (Color online) Stability of the spin-chiral phase and 
flatness of the topological bands depending on parameters of 
the Hamiltonian. In (a), shaded areas in the i^-AjT plane 
indicate a spin-chiral ground state Fig. |2^a,b) for U/t = 12, 
white areas have a different ground state. Shading indicates 
the figure of merit M for the flatness of the upper chiral sub- 
band, bright thick lines bound the region with M > 10. (b) 
shows M depending on tdd for selected sets of U and Ajt- 
Where the "Mott gap" , which separates the flat topologically 
non-trivial band from the upper Hubbard band, becomes very 
small, M is determined by the minimal gap separating the 
band of interest from other bands. J = U/A and t = 1 were 
used in all cases. 
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FIG. 4. (Color online) FQH state induced by NN Coulomb 
repulsion V in the effective one-band model Eq. [T] (a) Energy 
depending on total momentum k for several values of V/t. (b) 
Energy for V/t = 0.2 depending on a flux 4> y added whenever 
an electron goes once around the whole lattice in y direction. 
Each addition of (f> — 2n leads to an equivalent state, 6ir 
to the same state. The Chern numbers associated with the 
three low-energy states are almost exactly 2/3 for V/t = 0.2. 
Lattice size is 4 x 6 sites (12 two-site unit cells), parameters 
in Eq. [T] are t\ = 0.27t and £3 = — 0.06t, giving bands with 
M w 13 and a gap of 0.89*. The filling of the flat band is 2/3. 



Kondo-lattice model (KLM). Low-energy configurations 
minimize onsite interactions and thus contain two or 
three electrons per site, with parallel spins due to Hund's 
rule. In order to additionally minimize the crystal-field 
energy, the e' g levels will always be half filled and form 
an effective spin, while any holes will be found in the a± g 
sector. The electrons in the partially filled a\ g states can 
delocalize with an isotropic hopping t a i g = (2f + tdd)/3, 
however, their spin must remain parallel to the local e' g 
spin. In the low-energy limit, each site can thus be de- 
scribed as a spin coupled to a charge degree of freedom 
and we arrive at the situation described by the KLM in 
the limit of strong Hund's rule coupling. Our numeric 
mean- field results corroborate this picture, see Fig.[2|c), 
where the e' g levels are found far below the chemical po- 
tential. The KLM supports spin-chiral phases on many 
frustrated lattices like the triangular [3"51 I4TJ1 EPH 135] . py- 
rochlore [46], and face-centered cubic [47] lattices. 

In addition to processes within the low-energy Hilbert 
space, virtual excitations involving high-energy states 
can be taken into account in second-order perturbation 
theory. This leads to (i) effective longer-range hopping 
of the ai g electrons and (ii) an effective antiferromag- 
netic superexchange between the e' g spins. The latter 
stabilizes the spin-chiral pattern |45j and is due to ex- 
citations into the upper Hubbard/Kondo band. When 
it is suppressed for U > 24|i|, the ground state conse- 
quently becomes FM, as in the KLM with a large Kondo 
gap [ID] |33]- Nevertheless, the exotic spin-chiral state is 
remarkably stable in the present t2 g system considering 
its sensitivity to Hund's coupling in the KLM [40] . 

The effective longer-range hopping of a\ g elec- 



trons involves processes via excitations into the upper 
Kondo/Hubbard band (oc 1/J and oc l/U) as well as vir- 
tual excitations of e' g electrons into a\ g states (oc 1/Ajx), 
for details see EH Second-neighbor hopping oc 1/J does 
not significantly modify the low-energy bands and drops 
out completely in the limit of a large Mott/Hubbard gap, 
but third-neighbor hopping £3 is crucial in cancelling the 
dispersion coming from NN hopping t\ for one of the 
bands 35j . The simplest description of the effective low- 
energy bands around the Fermi level is thus 

ffefl(k) = 2ti ^ a 3 cos ka^ +2t 3 ^2 a ° cos 2ka j > C 1 ) 

where &j (j = 1, 2, 3) denote the unit vectors on the tri- 
angular lattice. Pauli matrices cr- 7 and unit matrix er° 
refer to the two sites of the electronic unit cell in the chi- 
ral phase (3S]. Formally, this describes electrons moving 
in a constant (and very strong) magnetic field with a flux 
of 7r/2 threading each triangle of the lattice [38] . 

FQH groundstates of an effective spinless one-band 
model. — We now address the impact of NN Coulomb in- 
teraction V J2u j) niUj on the fractionally filled flat band. 
The spin-chiral state can only be expected to remain sta- 
ble for densities close to 2.5 electrons per site, i.e., low 
doping factions v of the flat band [40 . FQH states corre- 
sponding to such low fillings arc generally separated from 
the rest of the spectrum by only a small gap, making their 
analysis on finite-size clusters difficult [48]. Here, we use 
Lanczos exact diagonalization (20] [35] |43l [48] to study 
a number of simple filling fractions (1/3, 2/3, 1/5, and 
2/5) available on accessible lattices and consistently find 
V to induce signatures of a FQH state. It is thus plau- 



4 



sible that the FQH behavior discussed next persists to 
fractional fillings in a low doping range of the spin-chiral 
state. 

As an example, we present here the case of 16 elec- 
trons on a 4 x 6-site cluster of the model Eq. ([I]), a fill- 
ing that would correspond to 2.6 in the original three- 
orbital model. After a particle-hole transformation, it 
corresponds to 2/3 filling of the nearly flat band. Fig- 
ure]^ a) shows that with increasing V, three low-energy 
states split off from the rest of the spectrum. Inserting 
a magnetic flux 4> y = 2ir interchanges the three states, 
4> y = 6tt recovers the original situation, see Fig. [Itb), 
as reported for other systems [23 GH 03]. The Chern 
number C is evaluated by integrating the flux-dependent 
Berry curvature Q n (<fi x , <j) y ) (obtained by the Kubo for- 
mula [351 ESI ISO]) over the square < <j> x ,<i>y < 67r. For 
V = 0.2|t|, the three low-energy states have Chern num- 
bers within 1% of the expected C — 2/3, a deviation well 
within the limits of reported finite-size effects [20] . 

Conclusions. — The possibility of a spontaneous FQH 
effect without a magnetic field is currently hotly dis- 
cussed, and various models have been suggested [TSl - 
[2171 |2"3TI2"6"1 14"3"] . However, an experimental realization ap- 
pears challenging, as the necessary topological character 
and the flatness of the bands need to be carefully engi- 
neered in previous proposals. We have shown here that 
bands with the desired properties emerge spontaneously 
for wide parameter ranges in strongly correlated ti g or- 
bitals on a triangular lattice, and that these bands sup- 
port FQH ground states. Both t^ g systems and triangu- 
lar lattices occur in various TM oxides, and signatures of 
the unconventional integer QH state have been reported 
for a triangular-lattice palladium-chromium oxide |51j . 
This harbors the prospect that a suitable material can 
be synthesized in this highly versatile material class. As 
such a material is by default strongly correlated, one also 
naturally expects an inter-site Coulomb repulsion that is 
strong enough to stabilize spontaneous FQH states in the 
absence of a magnetic field. 

This research was supported by DFG (Emmy-Noether 
program; SK and MD) and the Interphase Program of the 
Dutch Science Foundation NWO/FOM (JV and JvdB). 

The appendix contains Supplemental Material: 

One-particle terms of the multi-orbital t^g 
Hamiltonian 

The multi-orbital kinetic energy is 

H ^= E ^/4,a,.<w+ H - c ' ( 2 ) 

where c\ aa {c iaa ) creates (annihilates) an electron on 
site i, in orbital a and with spin a. (i,j) denotes nearest- 



neighbor (NN) bonds, a and (3 denote the orbital. Us- 
ing as basis states the xy, xz, and yz orbitals shown in 
Fig. 1(a) of the main text, the orbital- and direction- 
dependent hopping parameters t"'? are given by the ma- 
trices 

(t dd 0\ /0 A /0i 0\ 

fi = * , f a = t dd , f 3 = [t 
V t OJ \t 0/ \0 t dd J 

for NN bonds along the three directions ai, &2, &3 as 
illustrated in Fig. 1 of the main text. The transformation 
into the {a% g , e' + , e' _} can be found, e.g., in Ref. [33] . 

Mean-field approximation 

Onsite interaction is described by Kanamori param- 
eters U (U') for Coulomb repulsion between electrons 
in the same (different) orbitals as well as ferromagnetic 
Hund's-rule coupling between electrons in different or- 
bitals. The relation U' = U — 2J is used here, "pair- 
hopping" J' is left out, because it drops out of the mean 
field decoupling 

Hint ~ U ^ (( n i,a,t) n i,ca + n i,a,t ( n i,a,i)) 
i j a 

+ (U'- J/2) + ni, a (n i:0 )) 

— 2 J 22 (( m i,a) m ^+ m i,aK(3)) 
i,a<P 

i,a i,a</3 

+ 2J^2(rm, a )(m,p), (3) 

i,a</3 

where i labels the site, a and (3 orbitals. n^ a ^ — 
c lacr c iaa ^ s the density operator. We keep here only 
expectation values for diagonal operators, i.e., only 
(TH, a ,<r) = ( c l. a ,a c i,a.a) EZ1EH]> which reduces the spin 
to its z-component m^ a — (ni, Q ,t ~ n i,a,i)/2- In order 
to treat non-collinear spin patterns, one has to allow for 
a site-dependent spin-quantization axis given by angles 
0i and (f>i . The change in quantization axis from site to 
site manifests itself in a complex phase for the hopping 
terms, [39 which is between sites i and j 

Q°f = cos ^ cos ?l + sin ^ sin ?L e -^*i-*i) 
13 2 2 2 2 

Qt4 = cos °A sin ho-Hi - C os *l sin ^e"* (4) 
« 2 2 2 2 y ' 

where fiy (^Jy ) modulates the hopping of an elec- 
tron with spin parallel (antiparallel) to the chosen spin- 
quantization axis. In the site-dependent quantization, 
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spin is not conserved and there are spin-mixing hoppings 
with fijj given by the complex conjugate of ■ 

We use numerical optimization routines to find the spin 
pattern with the lowest energy among all orderings with 
unit cells of up to four sites, including all patterns consid- 
ered in Ref. [ID] of the main text. In each step, the mean- 
field energy is calculated self-consistently for a lattice of 
16 x 16 (four-site unit cell) or 24 x 16 (three-site unit cell). 
(For selected points in parameter space, we also used 
larger lattices and did not find a significant difference.) 
In order to minimize the impact of our approximations 
on the symmetries of the orbital degrees of freedom, we 
perform the mean-field decoupling in the {ai g , e' g+ , e' g _} 
basis, where the symmetry between the half-filled e' g + 
and the quarter-filled a\ g orbitals (for the fillings dis- 
cussed here) is already broken by the crystal field. We 
verified that decoupling directly in the {xy, xz, yz} basis, 
where all three orbitals have the same electronic density, 
leads to qualitatively identical and quantitatively very 
similar results. 



Effective one-band model and exact diagonalization 

The mapping to the effective one-band model is most 
easily carried out in the Kondo-lattice picture, where the 
localized spins are assumed to consist of the e' g elec- 
trons. Without a magnetic order, the a\ g orbital has 
an isotropic hopping t a i g — (2t + tdd)/3, i-e., the same 
along all three directions on the triangular lattice, but 
in the spin-chiral phase, this hopping is modulated by a 
direction-dependent Berry phase Eq. Q. The electronic 
unit cell of the spin-chiral pattern has two sites [35] , and 
the Berry phases can then be expressed in terms of Pauli 
matrices as given in the main text. The absolute value of 
the NN hopping is renormalized to t\ = (2i + tdd)/3V3. 

Corrections to this simplest approximation can be ob- 
tained by second-order perturbation theory, which yields 
longer-range hopping processes mediated by virtual exci- 
tations. We are first going to discuss processes within the 
lower Kondo/Hubbard band, where an electron from the 
e' g levels, which are half filled in the low-energy Hilbcrt 
space, is excited into an empty a\ g state in the virtual in- 
termediate state, involving an excitation energy Aj^. In 
a second step, an electron from a different occupied a\ g 
state can take the empty place in the e' g orbital, which 
corresponds to an effective hopping. There are two pos- 
sible hopping paths connecting pairs of either nearest or 
next-nearest neighbor (NNN) sites and it turns out that 
the corresponding effective hoppings drop out in the spin- 
chiral phase because the Berry phases for the two paths 
interfere destructively. Third-neighbor sites, on the other 
hand, are only connected by a single path and the com- 
bined Berry phase fijj^fij^ only renormalizes the effec- 
tive hopping by a factor of 3, because the spins at sites 
i and k are parallel in the spin-chiral phase. One thus 
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FIG. 5. (Color online) Comparison of the low-energy bands 
obtained in mean field theory for the spin-chiral state to those 
obtained in the strong-coupling Kondo-lattice model, includ- 
ing effective hoppings in second-order perturbation theory. 
Parameters in mean field are U/t = 12, J/t = 3, Aj-r/i = —6, 
and tdd = 0, as in Fig. 2 of the main text. The different 
dispersions in the upper band are due to the fact that the 
perturbation theory assumes that the second-order-hopping 
path is always available. As the ai g level is partly filled, some 
of these paths are blocked by the Pauli principle, this effect 
is included correctly in mean field. Taking this into account 
phenomenologically by reducing tz by 20% improves the fit. 
This was done for the curve referred to as u t3 only", where 
only nearest and third-neighbor hopping were included, while 
second-neighbor processes were dropped. 



obtains an effective hopping t 3 = —2(t — tdd) 2 / (27 Ajt), 
which is the only second-order correction for the limit of 
infinite Hund's rule. In this limit, £3 flattens the lower 
of the subbands for the fillings discussed here [53] . 

In the more realistic case of strong but finite Hund's 
rule coupling, there are additional processes where the 
virtual excitation involves an electron in the upper Hub- 
bard/Kondo band. The corresponding excitation en- 
ergies then contain Hund's coupling J and the effec- 
tive hoppings from site i to k involve the Berry phases 
^If^fk: see E< F For tne parameters of Fig. 2 of the 
main text, the chiral bands around the Fermi level are 
compared to this second-order treatment in Fig. [5] NNN 
hopping does here not drop out, and NN hopping is also 
slightly renormalized, however, both these processes have 
only a small impact on the low-energy bands. Again, we 
find third-neighbor hopping t$ to flatten the dispersion of 
one of the bands. Since excitations into the upper Hub- 
bard/Kondo band involve electrons rather than holes as 
before, the sign of the effective £3 is reversed, and it is 
the upper chiral subband that is flattened. 

In order to investigate the FQH groundstate, we used 
Lanczos exact diagonalization to study the Hamiltonian 
given by the kinetic energy of the effective one-band 
model, Eq. (1) of the main text, and NN Coulomb repul- 
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sion VJ2(ij) n i n j- Hopping parameters t\ = 0.27 and 
= —0.058 were used, giving a dispersion similar to 
Fig. [5} NN bonds are defined on the original tri- 

angular lattice and V acts both between the two sites 
within one unit cell and between NN sites belonging to 
different unit cells. As mentioned in the main text, the 
spin-chiral state can actually only be expected to be sta- 
ble for low doping of the flat band [40] , which is close to 
half filling for the effective one-band model. FQH states 
corresponding to such small fillings v tend to have smaller 
gaps than those for large v, and the low-energy manifold 
giving the quasi-degenerate FQH states contains more 
states. On the small clusters that we can study with ex- 
act diagonalization, eigenenergies always have spacings 
between them, as an illustration see the V = energies 
in Fig. 4(a) of the main text, which would form a con- 
tinuous band in the thermodynamic limit. It is thus far 
harder to reliably resolve a small gap than a larger one, 
and it is moreover highly desirable that we can study the 
system on at least two lattice sizes in order to see a gap. 
This severely restricts our access to very low dopings. We 
thus study several filling fractions corresponding to "sim- 
ple" FQH states. In all cases where we find a low-energy 
manifold to separate from the rest of the spectrum, the 
states of this low-energy manifold shows signatures of 



FQH behavior, which is thus a very robust feature of the 
doped flat band. 

Inserting a flux (<f> x ,(j) y ) means that electrons gain a 
phase e l ^ x (e 1 ^) for going once around the whole lattice 
in x-(y-) direction. This is implemented by changing the 
hopping ty from site i = i x &i + i y a.2 to site j = j x &i + 
jy& 2 to 

tij^iije 1 ^^^^), (5) 

leading to a flux-dependent Hamiltonian H c ^{<j) Xl <j) y ). In 
the case of v = 1/3, we find three low-energy states 
separated from the remaining spectrum by a gap as in 
the v = 2/3 case discussed in the main text, both on a 
6x6 and a 4 x 6 system. For 6x6 sites, however, all 
three low-energy states have total momentum (0, 0) for 
{(j) x ,(j)y) — (0,0). Due to finite-size effects, the states do 
then not cross upon flux insertion (35], but avoid cross- 
ings. For the smaller 4x6 system, the three low-energy 
states have different total momenta, and this good quan- 
tum number allows us to clearly resolve their crossing 
when we insert a flux <j) y , even on a finite system. 

The Chern numbers were evaluated by integrating the 
Berry curvature Q n ((f> x ,4> y ) over the square < <f> x , <p y < 
Qir. n n ((f) x ,<f)y) was obtained by the Kubo formula [13 

ED] 



n n ( 



iL x L y S ^ 



(n\ aH it M \ n ')( n '\ dH it:' M \K) 



1 9H a - 



n) 



(En - £n') 2 



(6) 



where n' and n label eigenstates with energies e n > /„ and 
dH e s{(j) x ,(j)y)/d(j) x i v are current operators. 
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